Distribution of Partial Neutron Widths for Nuclei Close to a Maximum of the 

Neutron Strength Function 



Hans A. Weidenmiiller 1 

1 Max-Planck-Institut fur Kernphysik 
Heidelberg, Germany 

For nuclei near a maximum of the neutron strength function, the secular dependence on energy E 
of s-wave partial neutron widths differs from the canonical form \[E~. We derive the universal form 
of that dependence and show that it is expected to significantly influence the analysis of neutron 
resonance data. 
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Purpose. The Porter-Thomas distribution [l[ is one of 
the key predictions of Random-Matrix Theory (RMT). 
The reduced partial neutron widths (simply "neutron 
widths" in the sequel) of compound-nucleus (CN) res- 
onances are predicted to follow a x 2 distribution with 
v = 1 degrees of freedom. That prediction was recently 
tested with unprecedented accuracy {2|. The authors 
scattered slow neutrons on several Pt isotopes, thereby 
measuring sequences of CN resonances over energy inter- 
vals of up to 20 keV length, and obtained sets of widths 
for s-wave neutrons containing up to 450 data points. 
Reduced neutron widths were obtained by rescaling the 
measured widths by the factor f 2 {E) = \f~E~ with E 
taken at resonance energy. That factor is supposed to 
take account of the secular variation with energy of the 
widths for s-wave neutrons. Using a cutoff procedure to 
minimize p-wave background and a maximum-likelihood 
analysis, the authors concluded that the validity of the 
Porter-Thomas distribution must be rejected with a sta- 
tistical significance of at least 99.997 per cent ] 2\. That 
result calls into question earlier successful tests of RMT 
in nuclei (for a review, see Ref. Q). More generally, it 
questions whether at excitation energies of several MeV 
nuclei are correctly described as basically chaotic sys- 
tems, a view widely held so far. Not surprisingly, the 
result announced by Koehler et al. has found wide at- 
tention [3]. 

For the isotopes of Pt investigated in Ref. Q , the neu- 
tron strength function (the ratio of the average neutron 
partial width and the mean resonance spacing) is strongly 
enhanced. The enhancement facilitates the separation 
of s-wave resonances from the contamination of small 
p-wave resonances. To the best of our knowledge, the 
theoretical implications of that enhancement for the sec- 
ular dependence of neutron widths on energy, i.e., for the 
function f 2 {E) introduced above, have never been inves- 
tigated. In this Letter we show that the enhancement of 
the strength function implies that f 2 (E) carries an addi- 
tional energy dependence on top of the V^E-dependence 
mentioned above. That dependence has never been taken 
into account in any analysis of data on neutron widths. 
It is not clear whether the conclusions drawn in Ref. |2j 



will be upheld when that dependence is accounted for. 

Maxima of the neutron strength function are due to 
a resonance close to threshold (more precisely: to a vir- 
tual state as defined below) or to a weakly bound state 
of the nuclear single-particle potential for s-wave neu- 
trons. Such maxima occur systematically with increas- 
ing mass number (i.e., increasing radius of the nuclear 
single-particle potential) when a bound s-wave state is 
being formed at threshold. (For the Pt isotopes that is 
the 4s-state in spectroscopic notation). The enhance- 
ment of neutron widths due to the nascent bound single- 
particle state is strongest for CN resonances close to neu- 
tron threshold and, being a threshold effect, dies out with 
increasing separation of the CN resonances from neutron 
threshold. Such weakening of the enhancement causes 
the additional energy dependence of f 2 {E) and forms 
the topic of the present paper. 

To make the paper self-contained we begin with a brief 
account of single-channel CN scattering theory. We show 
that f(E) is determined by the projection of the neutron 
single-particle scattering wave function onto the nuclear 
volume. We calculate that projected part of the wave 
function for an attractive square-well potential and de- 
termine f(E) for both, a single-particle resonance close 
to threshold and a weakly bound state. We show that 
the enhancement is a generic feature and occurs like- 
wise in more realistic nuclear single-particle potentials 
that differ from a square well. Whenever the s-wave 
strength function shows significant enhancement, f 2 (E) 
differs significantly from y/E even over energy intervals 
as small as 20 keV (the interval used in Ref. |2j). 

Single-Channel CN Scattering. We deal with a single 
channel (the s-wave neutron channel). The energy E is 
positive and E = denotes the threshold. In the shell- 
model approach to nuclear reactions [H , the channel wave 
function \e is the antisymmetrized product of the wave 
function of the target nucleus in its ground state and the 
neutron single-particle scattering wave function ^ e (f) ■ 
The latter depends only on the radial coordinate; r and is 
chosen real. The numerous many-body resonances typ- 
ical for CN scattering are due to N 3> 1 orthonormal 
quasibound states where /i = 1, 2, . . . , N. These in- 
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teract with each other through the ./V-dimensional real 
and symmetric Hamiltonian matrix (fi\H\i>) = 7J M „ and 
are coupled to the neutron channel by real Hamiltonian 
matrix elements 

W^E) = ( X e\H\h) . (1) 

The unitary scattering amplitude S(E) that takes into 
account the presence of the N CN resonances has the 
form d i] 

S(E)=exp[2i5}(l-2iK ]T W p (E)D^(E)W l ,(E)\ . 

(2) 

Here 6 is the s-wave potential scattering phase shift, and 
D lu ,(E) = E8 lu ,-H llv -F lu ,(E) + i*W M (E!)W v (E) (3) 
where 

r°° t> 
F^(E) = J dE' ——W^E')W„{E') (4) 

and where V denotes the principal- value integral. The 
last term in Eq. ([3]) is proportional to the width matrix 
and descibes the instability of the CN resonances due 
to their coupling to the channel. The matrix F^{E) in 
Eq. (fj| accounts for the shift of the resonances due to that 
coupling. Eqs. J5J to JD provide a very general frame- 
work for s-wave neutron scattering in the presence of CN 
resonances. The assumption that the CN resonances are 
described by RMT is implemented by postulating that 
the matrix H pll is a member of the Gaussian Orthogonal 
Ensemble (GOE) of random matrices. 

Near neutron threshold the CN resonances are iso- 
lated, and we use the diagonal representation H^ v = 
J2 P O pp E p O pv where M „ is orthogonal and where E p are 
the eigenvalues of H^. We dchne W^E) =Y^ V ®iivW v 
and F by Eq. ((4]) with all Ws replaced by Ws. For iso- 
lated resonances the coupling between the eigenvalues E p 
is negligible (both F M „ and W^Wv are taken to be diag- 
onal), and S(E) takes the form 

S(E) = exp[2i6] (l - 2in W^ (E) (E - S^W^E)) 

(5) 

where the complex resonance energies £ ^ are given by 

£» = E p + F^ - inW 2 . (6) 

The neutron partial width amplitude of the res- 
onance located at £ M is given by V2ttW p (E) = 
V%xYlv @ttvW p (E). In order to remove any secular en- 
ergy dependence from the matrix elements Wp(E) we 
write 

W^E) = f(E)V M . (7) 



By construction, the amplitudes do not carry any 
secular dependence on energy E. GOE predicts that in 
the limit of infinite matrix dimension, the elements Op U 
and, therefore, the amplitudes J2 V O pv V v are Gaussian- 
distributed random variables. As a consequence, the re- 
duced neutron widths 2ttW 2 (E)/ f 2 (E) are predicted to 
follow the Porter-Thomas distribution. To test that pre- 
diction, we must determine the function f(E) for nuclei 
in the vicinity of a single-particle s-wave resonance close 
to threshold and for a weakly bound single-particle s- 
wave state. 

Poles of the Single-Particle Scattering Amplitude. 
Eq. (TT|) shows that the entire energy dependence of 
Wp(E) is due to xe- We recall that xe is the anti- 
symmetrized product of the ground-state wave function 
of the target nucleus and the real s-wave scattering wave 
function ^£;(r). For the matrix element W P (E), only 
the projection of ^e onto the nuclear volume is relevant. 
Thus, the energy dependence of Wp(E) is determined by 
the enery dependence of ^^(r) for r < R where R is the 
nuclear radius. 

The scale in energy over which the radial dependence 
of a single-particle scattering wave function changes, 
is given by the typical distance between bound s-wave 
single-particle states. In a heavy nucleus, that distance 
is about 10 MeV and, thus, very large compared to the 
typical energy scale over which resonance data are taken. 
(InRef. 0, that scale was 20 keV). Therefore, the radial 
form of the neutron s-wave scattering function in the nu- 
clear volume r < R can safely be taken as independent of 
energy, and we focus attention on the energy-dependent 
amplitude f(E) of that function. In the vicinity of the 
threshold energy (E = 0), f(E) is strongly enhanced 
whenever a single-particle s-wave resonance or a bound 
s-wave single-particle state occurs close to threshold. 
Both resonance and bound state manifest themselves as 
poles of the unitary single-particle potential scattering 
amplitude s(E) — exp[2iS] in Eq. ([2|). To understand 
qualitatively what happens we recall some properties of 
the poles of s(E). The potential scattering wave function 
^e(t) at energy E depends asymptotically on wave num- 
ber k where h 2 k 2 /(2m) = E, and it is useful to consider 
the distribution of poles of the scattering amplitude s(E) 
in the complex /c-plane (rather than the complex energy 
plane). We accordingly replace s(E) by s(fc). 

For a square- well potential, poles of s(k) have been 
studied in detail 0,0. Poles of s(k) occur either pair- 
wise or as single poles. Pairs of poles lie below the 
real fc-axis and occur symmetrically to the imaginary 
fc-axis. Such pairs are, thus, located at fco and at — fcg 
with ^s(ko) < —i/a (where a is the radius of the square- 
well potential) and Ji(fco) > 0. Single poles lie on the 
imaginary fc-axis. Poles on the positive imaginary fc-axis 
correspond to bound states while poles on the negative 
imaginary fc-axis are referred to as virtual states. We 
visualize the motion of the poles in the complex fc-plane 



3 



as the depth Vq of the square-well potential is increased. 
For a very shallow potential, there are no bound states 
and no poles on the imaginary fc-axis. All poles lie far 
below the real /c-axis, and significant resonant behavior 
is absent. As Vq is increased, the first pair of poles ap- 
proaches the negative imaginary axis from opposite sides. 
The two poles coalesce at k = —i/a. Then one pole moves 
up and the other pole moves down the imaginary /c-axis. 
Significant resonance behavior of s(fc) is caused only by 
the upward-moving pole, first as a virtual state and later 
as a weakly bound state. (For a = 6 fm, a realistic value 
for the radius of a heavy nucleus, the point k = —i/a has 
a distance in energy of about 0.5 MeV from threshold). 
As Vo is increased further, resonance enhancement sub- 
sides. The pattern repeats itself as the next pair of poles 
approaches the point k = —i/a and the 2s-state is pulled 
into the potential, and so on. 

Enhancement factor for the square-well potential. We 
work out the resulting enhancement of the amplitude 
f{E) for the 4s-state although the analysis and result 
are exactly the same for any s-state near threshold. The 
real scattering function >f e(t") is normalized to a delta 
function in energy. An elementary calculation shows 
that for r < a we have <J e(v) — f(E)sm(nr)/r where 
h 2 n 2 /(2m) = E + V and where 



f(E) 



2(ka) 



wkh' 



J (ka) 2 sm 2 (ka) + (na) 2 cos 2 (na) 



(8) 



As remarked above, the function sin(Kr) changes very 
slowly with energy, and attention is focused on f(E). 

A weakly bound 4s-state with energy binding Eq < 
exists if the condition [(fca)/(fea)] tan(fca) = —1 is met 
for (na) = (7n/2) + e and e < 1. Then E = 
~(77r/2) 2 (h 2 e 2 )/(2ma 2 ). Expanding (na), sin(Ka) and 
cos(na) in Taylor series around na = 7n/2 and keeping 
only lowest-order terms, we find 



E+\E \ 



(9) 



We have suppressed energy-independent factors. The 
factor \f~E~ is universal for s-wave scattering near 
threshold. The factor (E + \EqD~ 1 describes the en- 
hancement due to the weakly bound single-particle s— 
wave state. For the virtual state the condition reads 
[(fea)/(Ka)]tan(Ko) = +1. That yields for f 2 (E) the 
same form as in Eq. (J9j) except that now Eq is the energy 
associated with the virtual state on the negative imagi- 
nary /c-axis. 

We have expanded {na) in powers of (ka) around the 
location of the pole of s(E), and we have kept only 
the term of zeroth order. The term of next order is 
(ka) 2 /(7n). For < E < 100 keV and a = 6 fm, that 
term is not larger than 0.01, and we expect the pole ap- 
proximation in Eq. Q to be excellent. 



Discussion. Although derived specifically for a square- 
well potential, the factor given in Eq. ((9]) is universal and 
describes for r < R the enhancement of the square of the 
s-wave single-particle scattering wave function also for 
other, more realistic forms of the single-particle poten- 
tial. Indeed, the pattern of movement of the poles of s(k) 
in the complex /c-plane versus depth of the potential is 
the same for all potentials that lack a barrier. Pairs of 
s-wave poles occur some distance below the real /c-axis 
and symmetrically to the imaginary axis and cannot give 
rise to significant resonance enhancement. Therefore, the 
value of k where such pairs coalesce on the negative imag- 
inary fc-axis, although different from that of the square- 
well potential, amounts to at least several 100 keV. Signif- 
icant resonance behavior is again due to the pole of s(k) 
that moves up on the imaginary fc-axis from the point 
of coalescence, first as a virtual and then as a weakly 
bound state. Potential scattering theory [j| shows that 
at a pole of s(k), \t_E(r) is singular. For k positive and 
close to a pole on the imaginary /c-axis, f 2 (E) has the 
form of Eq. ©. 

The universality of our result is displayed by the fact 
that f 2 (E) in Eq. © depends only on the energy of the 
virtual or weakly bound state and is independent of any 
detailed properties of the potential. (The value of Eq in 
Eq. ((9]) does, of course, depend on the potential). The 
denominator in Eq. ((9]) is obviously a special case of the 
universal Lorentzian factor [(E — Eq) 2 + (l/4)r 2 ] _1/ ' 2 de- 
scribing resonance enhancement and applies when the 
resonance is located below threshold so that E is nega- 
tive and r vanishes. We also note that the enhancement 
factor in Eq. ((9]) is similar to that due to a doorway state. 

We conclude that for a virtual or a weakly bound neu- 
tron s-wave state, the partial neutron widths carry the 
universal enhancement factor given in Eq. The en- 
hancement factor affects both, the neutron strength func- 
tion and the secular energy dependence of the widths of 
neutron s-wave resonances. To see what happens to the 
neutron strength function, we keep E (or k) fixed, posi- 
tive, and slightly above threshold, and we increase mass 
number A, thereby increasing the radius of the poten- 
tial. A pair of s-wave neutron resonances approaches 
the imaginary /c-axis, eventually giving rise to a virtual 
and, later, to a bound state. For the virtual state, the en- 
ergy | S 1 decreases monotonically toward zero and then, 
as the virtual state turns into a bound state, increases 
monotonically from zero. The result is first a steady 
increase of the enhancement factor (|9|) and then, after 
the virtual state becomes a bound state, a decrease of 
that factor. Taken together, that causes the characteris- 
tic maximum of the strength function. Conversely, when 
neutron widths are measured near threshold for a se- 
quence of s-wave resonances at some fixed value of A 
where the strength function is enhanced and, therefore, 
the enhancement factor in Eq. §§§ is operative, the re- 
sulting secular energy dependence of the neutron widths 
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goes beyond the standard V^-dependence. 

Given the universality of the enhancement factor in 
Eq. (|9]), we ask how that factor is expected to affect the 
analysis of neutron resonance data. We take the work of 
Ref. [U as an example. We assume that the distance 
of \Eq\ from threshold is large compared to the aver- 
age resonance spacing. We recall that the data typically 
range over an energy interval of 20 keV. To predict the 
influence of the enhancement factor in Eq. (|9]) on the 
data, one would have to determine Eq with an accuracy 
of about 100 keV. That requirement precludes a theoreti- 
cal prediction using the neutron single-particle potential 
in heavy nuclei. That potential is not known to such 
an accuracy. Therefore, we use another estimate. The 
mass dependence of single-particle levels is governed by 
the factor A -2 / 3 . For a potential depth of several ten 
MeV and A as 200, a weakly bound single-particle level 
changes by about 150 keV when A changes by one unit. 
Therefore, we expect that for nuclei near the maximum 
of the neutron strength function, \Eq\ is of the order of 
100 keV. That is also consistent with the fact that for the 
square-well potential, the point of coalscence has a dis- 
tance of about 500 keV from threshold. Taking in Eq. © 
lE'ol = 50 keV (10 keV) as examples, we see that the res- 
onance enhancement factor (the denominator of expres- 
sion ©) changes by a factor 1.4 (a factor 3, respectively) 
over an interval of 20 keV starting from threshold. Thus, 
the resonance enhancement in Eq. @ may significantly 
influence the analysis of data on neutron widths. 

In such an analysis, a first estimate of Eq can be ob- 
tained from Eq. ^ and the measured enhancement of the 
strength function. One may then consider E and the av- 
erage value of the reduced neutron partial widths as free 
parameters, and search for a best fit to a \ 2 distribution 
with v degrees of freedom for the reduced partial widths. 
Alternatively, it is possible to determine the enhancement 
factor directly, i.e., without using the pole approximation 
of Eq. (j9}. That can be done by solving numerically the 
radial Schrodinger equation for s-wave neutrons for a 
potential that is realistic in form and that possesses a 
virtual or a bound state close to threshold. The factor 
f(E) is then determined by the delta-function normaliza- 
tion condition for the solution of the radial equation. We 
stress that the single-particle nuclear potential to be used 
is not the optical-model potential (which describes the 
average neutron scattering amplitude and, thus, incorpo- 
rates the effect of the CN resonances) nor is it the real 
part of the optical-model potential (because the imag- 
inary part contributes to the real part via a dispersion 
relation). Rather, it is the pure shell-model potential for 
neutrons. Conclusions about the failure of RMT can be 
drawn only if such an approach yields a value for v that 
is significantly different from unity. Conversely, if agree- 
ment is obtained for v fts 1, that should make it possible 
to determine the energy Eq fairly precisely (with an error 
of perhaps not more than 100 keV) and, from there, the 



shell-model potential for neutrons in heavy nuclei with 
great accuracy. 

We mention in passing that the resonance enhance- 
ment in Eq. ([9]) also affects the shift matrix F llu (E) in 
Eq. (|4|. With the same approximations as used above, 
the level shift F MM in Eq. ((6]) is given by the product of 
the reduced partial width for the resonance at £ M and 
the factor J °° dE' Vf 2 {E')/{E^ - E'). The enhance- 
ment of the shift due to the pole of f 2 (E) in Eq. (|9]) may 
be considerable and have important implications for the 
analysis of spectral fluctuations. 

So far we have considered the case where \Eq\ is much 
larger than the mean spacing of the resonances. That is 
probably the typical case. If that condition fails (i.e., if 
l-E-ol < 5 keV or so), it is not justified to consider the 
neutron widths as energy-independent constants. Then 
it is essential to display the full energy dependence of 
S(E) explicitly. We use Eqs. (J7|) and (j9|) and rewrite the 
last term in Eq. ([5]) as 

-2wr^ v^V^[(E+\E \)(E-E^-F^)+in^EV 2 ]-% . 

(10) 

Here = X^^W^- Eq. (ITU|) is qualitatively different 
from Eq. ([3]) with energy-independent widths. The for- 
mulas of i?-matrix theory commonly used for the analysis 
of neutron cross-section data bear a close analogy to the 
latter equation [5|. Eq. (fTO)) shows that these formulas 
cannot be used when E is very close to threshold. 

In summary, we have shown that a substantial en- 
hancement of the neutron strength function necessarily 
implies a significant energy dependence of the neutron 
partial widths for resonances in the vicinity of neutron 
threshold. We have derived the universal form of that en- 
ergy dependence. Conclusions about the validity of RMT 
predictions can reliably be drawn only when that depen- 
dence is taken into account in the analysis of neutron 
resonance data. 
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